Parton Energy Loss at Twist-Six in Deeply Inelastic e-A Scattering 
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Within the framework of the generalized factorization in pQCD, we investigate the multiple par- 
ton scattering and induced parton energy loss at twist-6 in deeply inelastic e-A scattering with the 
helicity amplitude approximation. It is shown that twist-6 processes will give rise to additional nu- 
clear size dependence of the parton energy loss due to LPM interference effect while its contribution 
is power suppressed. 



Jet quenching, or parton energy loss induced by multi- 
ple scattering in high-energy nuclear collisions has been 
proposed as a very sensitive probe of the properties of 
the hot and dense medium P, which recently has 
given a compelling theoretical explanation for many ex- 
citing experimental phenomena observed at RHIC, such 
as the suppression of high px single hadron produc- 
tion, and the disappearance of away-side two hadron 
correlation j^. To investigate the radiative energy loss, 
many theoretical methods have been developed so far 
|llIlllllllilISII3- Recently, based on the general- 
ized factorization theorem the twist expansion approach 
has been proposed to derive parton energy loss in terms 
of the modified parton fragmentation functions in nuclei 
0, 0, 01 , and provides the first evidence that the A^^^ 
dependence of the jet energy loss describes very well the 
HERMES data in e-A deeply inelastic scattering (DIS) 
|l2j | . This twist expansion approach has also been applied 
to discussing nuclear enhanced effects in the Drell-Yan 
process , heavy quark energy loss induced by gluon 
radiation in nuclei 0| and other processes |0, 0| . 

In this letter, we will extend the twist expansion ap- 
proach to study the parton multiple scattering at twist-6 
for e-A DIS in nuclei. In the framework of generalized fac- 
torization of high twist 0, ^3 , we will derive the semi- 
inclusive hadronic tensor at twist-6 with the approach of 
helicity amplitude approximation(HAA) and the cor- 
responding parton energy loss. Also we will show the nu- 
clear size Ra dependence of parton energy loss at twist-6 
due to the Landau-Pomeranchuk-Migdal(LPM) interfer- 
ence effect 19]. 

We consider the semi-inclusive process in e-A DIS, 
e(Li) + A{p) — > e(L2) + h{lh) + X. The differential 
cross section for this process has a form as 
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where p ~ [p+,0, Oj^] is the momentum per nuclcon 
in the nucleus with the atomic number A, q ~ L2 — 
Li — [—Q'^/2q~,q~,0±] is the momentum transfer, the 
Bjorken variable is defined a,s xb = /2p'^q^ , s = 
(p + Li)"^ and a^M is the electromagnetic (EM) couphng 
constant. The leptonic tensor is L^^ = 1/2 Tr(7 • £17^7 • 
L2lv) while the semi-inclusive hadronic tensor is defined 
as 



dHu 2 

X 2TTS'^{q + p-px - ih), 



X 
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here runs over all possible final states and = 
J2q ^q''i'qltJ.'4'q IS the hadronic EM current. 

In the parton model with coUinear factorization ap- 
proximation, the leading-twist contribution at the lowest 
order of single scattering can be factorized as, 



dW^ 

dzh 



J2 J dxff{x,^,j) 



Hj^J {x,p, 9) = Y Tr(7 • p7^7 • {q + xp)-/^) 
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where the momentum fraction carried by the hadron 
is defined as Zh = t'^jq^. Hjj^J is the hard par- 
tonic tensor, /i| and are the factorization scales for 
the initial quark distributions fq{x,^'\) in a nucleus 
and the fragmentation functions Dq^h(zh, ^i^), respec- 
tively. Including all leading log radiative corrections, the 
renormalized quark fragmentation function Dq^h{zhi M^) 
satisfies the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution equations pol |. 

In the nuclei the propagating quark in DIS will ex- 
perience additional scatterings with other partons from 
the nucleus. The rescatterings may induce additional 
gluon radiation which will cause the energy loss of the 
leading quark. Such phenomena are associated with the 
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FIG. 1: Sample diagrams for rescattering with gluon radiation 
in deeply inelastic e-A scattering, x, a;i, X2, x^, X4 and X5 are 
the momentum fractions carried by initial parton. 0, y, yi, 
2/2, ys and 2/4 represent the positions of the rescattering. 



high twist effect. Previous results H ES [H El El El 

have discussed the double scattering which gives twist-4 
contribution. Here we would like to go beyond double 
scattering and extend the study to investigate the twist- 
6 (triple scattering) processes as illustrated in Fig.Q] and 
their contribution to the parton energy loss. 

According to the generalized factorization theorem of 
high twist processes 0, 0, Esj, the leading twist-6 
(triple scattering) contribution with nuclear medium ef- 
fect to the semi-inclusive hadronic tensor can be formu- 
lated as 



"1 dz 



X / ^ dy^ dy^ dy^ dy^ Kgg^y ^Vi ^V2 ^V:i ^Va 
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4! dk;^d4 dkt^dk- 



HaviV ,2/1 ,^2 '^4 'fcT,P,g, Z) 



(4) 



where kx is the initial gluon's intrinsic transverse mo- 
mentum and T^^g is related to the parton correlation 
function of one quark and two gluons in a nucleus, 
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T'qgg[y ,2/1 ,^2 .2/4 ) = ^ (^1^9(0)7^ 

F+ [y^ ) F+ {y^)F+- [y^ ) (yf )^,{y-) \ A) .(5) 

Because of the collinear approximation, the tensor struc- 
ture of the triple scattering is generally the same as in 



the leading order single scattering, H^^, in Eq. Q can be 
expressed as 

H^Av^ ,Vi ,y2 ^Va ^kT,p,q,z)^ j dxHj^J {x,p,q) 



J" _____ 

x^^ (y ,2/1 ,2/2 '2/4 >^T,2:,p,9,z) . (6) 

At twist-4 H El El El El there are 23 cut-diagrams, 
while at twist-6 there are 201 cut-diagrams. To sim- 
plify our calculations, we adopt the helicity amplitude 
approximation with the soft gluon radiation Zg — 

(1 — 0) — > 0, where Zg is the momentum fraction carried 
by the radiated gluon and z is the momentum fraction 
carried by the leading quark. With helicity amplitude 
approximation, it is similar to the twist-4 processes that 
the leading order contribution comes from 49 central-cut 
diagrams for the twist-6 processes. 

For the contributions of the diagrams in Fig. ^ the 
results can be expressed as 

K) « ^(l-e'-^'"^0 

X pi{xB+XL)p'''y^+ixDP''' ivi +y3 -Vi) 
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JXLP'^y^ 



X gi{xB+XL)p'''y +ixDP'^(.yi -^2 +^3 -2/4 ) 

It is similar to the results of Eqs. iQ and ©, one can 
check that the contributions from 40 diagrams among the 
49 central-cut diagrams have the form 



— T 

H oc 



rjkj 



(9) 



where we introduced a symbol rj — X and Y repre- 
sent the longitudinal momentum fraction and the spatial 

coordinates, respectively. In Eq. Q, the form of Hf^^-^ 

corresponds to the case for = and i/^^^ in Eq. is 
the case for -q ~ \. Other 38 diagrams have the similar 
results. According to the above factorization formula, 
we can prove that the fourth derivative of the above ex- 
pression with respect to k^ vanishes at fcT=0 when we 
keep only the leading terms as It — > 0. Therefore, the 
contributions from these 40 central-cut diagrams to the 
hadronic tensor will vanish. As a result, only 9 central- 
cut diagrams as shown in Fig. |21 give non-zero contribu- 
tion to the hadronic tensor with helicity amplitude ap- 
proximation. 

For the central-cut diagrams in Fig. |2 the contribu- 
tions have the following form 



as 1 + 2% ^ 2 2 

— T. ; iDTT q;„ 

{It - krY 27r 1 - z 



^ ^i{x+XL)p*y +ixop^{yi -2/2 +^3 "Vi '>6{—y^) 
X %4"-2/3")^(2/r- 2/2^)^(2/" -2/rK- (10) 
For each diagram in Fig. [21 li can be expressed as 



(^g-H^£)-^-D)P+(j/4 -^3 )-iXLP'^y^ 
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FIG. 2; Diagrams of the twist-6 processes with gluon radiation and non-zero contribution to the hadronic tensor in deeply 
inelastic e-A scattering. 



_ gi(l-2:/(l-z))a:DP^(yr-yir)-*^-£-P'''f4 ) 

72 = g-«(l-z/(l-z))a^DP*(y3^-y4 ) 

X (^e-^^i-P^yi — g«^-DP+j/7/(i-z)-) 

_ g-i(l-z/(l-2))2:DP+(t/7-J/2")-«^i.P^(j'"-ar) 

_ g-«(l-2:/(l-2))2:DP+(yr^2'2") 

— g*(l~^/(l~^))=^DP^(!^4^"J'3^)~''^-C'P^J'4 ) 

— (^g-«2;i,p+a4" _ g'^-DP^ 1/4^/(1-2)) 

X (^e^(^"D^^ o)p'^ ivi LP''' {y^ -y2 ) 



, (11) 



) ,(12) 



(13) 

) ,(14) 
, (15) 



Jg = (^g-i(x"o-^D)p*{yi-y3)-ixLP*yi; 

_ g«(i-2/(i-2))2:DP+(y4"-y;7)-«a:i,p+i/;;; -j 

_ g-i(l-z/(l-z))a,£jp+(a7-y2")-"iP^(2'"~yr)) ^(16) 

= g-«(i-2/(i-2))2:DP+(ar-y2"+a3"-?'4r) 

X (^e'-^^LP'^Vi — g^^op^Vi /i'^^^-)'^ 

X (e 

Is = e 

X (g-':riP+(a--yr) _ g 

7g — g— '(l-2/(l-z))a:i3P+(y3"-y4 ) 

X {^g-^^LP^Vi — g^^-DP^J'4 /(l-z)) 

X (e 



i(l-2/(l-2;))2:£5p+(t;j"-j/2") 

!X£5p+J/4 /(I 



-ixip+(a -1/1 ) _ ^ixDP^iy )/(l-z)^ ^ ^29) 



The variables x^, and xu in the above expressions 
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are defined as 



XL 



2p+q-z{l - z) ' 



2p+q- 



xd 



2kT ■ £i 



2p+q^z 



(20) 
(21) 
(22) 



where £t is the transverse momentum of the radiated 
gluon, and z — (q/q~ is the momentum fraction car- 
ried by the final quark. Four terms of each expression in 
Eqs. (|11I19|) represent the contributions from different 
processes and their interferences. In the collinear Umit 
as ^ and fcy ^ 0, there will be a clear cancelation 
of these contributions due to the LPM effect. We will see 
later that the LPM effect will give an additional nuclear 
size dependence to the parton energy loss. 

Adding up all contributions together, we can derive 
the semi-inclusive hadronic tensor for triple scattering 
processes at twist-6, 



1 [ ^I67r2«2[i_ 



1 + ^ 



■A- 



+ <5(2-1)AT4(x,4)], 
where the " -t- " functions have a form 
Fiz) _ 



(23) 



dz- 

(1-^)h 



\^m^, (24) 
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F(z) being any function which is sufficiently smooth at 
z = 1. Here we also include the virtual corrections in 
Eq. (|23|1 . which is obtained by using the unitarity re- 
uirement and ensures the final result to be infrared safe 



The parton correlation function {x,xl) in Eq. H23() 
is expressed as 



rpA 



dy 
2-K 



dy^ dy^ dy^ dy^ e 



(25) 



x^(yi - y2 Wiy - Vi Wi-y^ Wiy4 ~ % ) 

xIixL)T'gfg{y-,y^,y2,y3,yi) , 
where function I{xl) is defined as 

i=l,2 j"=3,4 

X (1 -e-"^f^(^"-^-")). (26) 

Here Cij are associated with the color factors of each cut- 
diagram, and we have C14 — ||, C24 = C13 ~ —-§1 and 



C23 = j2' principle, the parton correlation function 



{x,xl) is not calculable and can only be measured 



in the experiments. However, under some assumptions, 
one can relate this parton correlation function to parton 
distributions and estimate its value [13, El The 
definition of AT^gg{x,e'^) in Eq. ^ is 



dz- 



[2Tf(x,XL) 



- {l + z^)T^(x,XL 



(27) 



Considering the high- twist contributions up to twist-6, 
we can express the hadronic tensor as 



dW^ 

dzh 



I dxf^{x,fij)Hj^Jix,p,q) 



X Dg^ft(Z/i,^^), 



(28) 



where the quark distribution function f^{x,^^) include 
also the twist-6 contributions. Up to twist-6 the modified 
quark fragmentation function Dq^h{zh, M^) has the form 



Dq^h{zh,^i^) = Dq^h{zh,^l^) + d''*^1^* ^\zh,^l'^) 

+D\__i, '{zh,p) 

= Dg^h{zh, ) + D' '{Zh,fl)+ 

Jo *T 



1 dz 



(29) 



" 27r ' - '~ 
+A7(*-Jf (z, X, XL,er)Dg_,n{znlz)] 



where quark fragmentation function D^^^^* ^\zh,^i'^) 
together with corresponding modified splitting functions 

^'lq-^qg~'^\z,X,XL,(T) a-Ild A'J^*^q~'^\z,X,XL,iT) 

twist-4 have been derived in Ref. H 0| , the third term 
(zft, /i^) in Eq. (|29|l is the additional corrections 
to quark fragmentation function at twist-6. In Eq. I|29|) 
the modified splitting functions at twist-6 are 



2(4^3 



3(4 + (fc2))2/,^(x,Mf) 



1+z' 

\i-zy 



-T''(x,XL)+S{l-z)AT^(xj'r)] , (30) 



A-f('^lf-'^\z,X,XL,fT) 



A7, 



(twist— 6) /'I 
q^qg \^ 



(31) 



Here, we replace l/i^ in Eq.^ with l/£'^{£'^ + (k^))^ in 
order to restore the collinear structure of the scattering 
amplitude [lo| |. 

In Ea. (|23(l for the hadronic tensor at twist-6, there is an 
additional factor as compared to the case at twist-4 
For large LPM effect becomes unimportant 
and one can neglect the interference contribution. In 
these processes, large final transverse momentum £^ ^ 
will lead to a suppression as compared to the 

case at twist-4 ^3,^3- On the other hand, when i'^ takes 
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an intermediate value, we can not ignore LPM effect. 
One can check in Eg. 1)26(1 . if — > 0, there will be a clear 
cancelation by the interferences. The LPM interference 
effect restricts the radiated gluon to have a minimum 
transverse momentum ~ Q"^ jA^I^ JJij. Therefore, the 
additional factor 1 jGp in Eq. H2^^|l as compared to the case 
at twist-4 has a suppression of Ra/Q^ for intermediate 
at twist-6 due to the LPM interference effect, here 
Ra ~ 1.12^^/'^ is the radius of the nucleus. 

As shown in Ref. [13,113, the fractional energy loss of 
the leading quark can be given by the fractional energy 
carried away by the induced gluon, the fractional energy 
loss of the leading quark at twist-4 can be expressed as 



, a,, 
dz — z 
2tt 



1 l+{l-zf alT^g{xB,XL) 



dz 




IHxb) 



,(32) 



here {xb ) is the quark distribution function up to the 
twist-4. For the twist-4 processes, the parton correlation 
fimction T^g{x,XL) contributes a Ra dependence [Tol. IT^ . 
After taking into account 1/^^ ~ Ra/Q^ resulting from 
the LPM effect, the nuclear correction to the fragmenta- 
tion function and the parton energy loss due to double 
parton scattering (twist-4) processes will then be in the 
order of alRA/ir ~ alR^/Q^. As shown in Ref . [l3,[l3, 
the order of the fractional energy loss for leading quark 
has a form 



XB 



(33) 



where xa = ^/tunRa with m^v is the nucleon's mass. 

At twist-6 it is the same as in the case of the twist-4 
the fractional energy loss of the leading quark is given by 
the fractional energy carried away by the induced gluon, 
then the fractional energy loss of the leading quark has 
a form 



2 I dz — z 



xA^j;i^lf-^Hz,XB,XL,fT) 

l + {l-z)^ 



del 



dz- 



l6TTalT^gg{xB,XL) 
3/^(xb) ' 



(34) 



here quark distribution function /^{xb) contains the 
contribution of the twist-6 processes. The parton cor- 
relation function T^g{x, x^) at twist-6 will contribute 
R\ dependence because T^g{x, x^) has an additional Ra 
dependence due to the extra gluon pair as compared to 
the case at twist-4. Here we adopt also the assumption 



[Tol IT^ IT^ that the positions of the two extra gluon 
field strengths are confined within one nucleon. As has 
been shown Ra/Q^ due to the LPM effect, com- 

bine R\ dependence of the parton correlation function 
T^g{x,XL), the nuclear correction to the fragmentation 
function and the parton energy loss due to triple parton 
scattering (twist-6) processes will then be in the order 
of a^R\/i'^ ~ a^R\/Q'^. The order of the fractional 
energy loss for leading quark can be expressed as 



g /twist — 6 



(35) 



One can see {Azg)twist-6 has asR\/Q^ suppression as 
compared to {A.Zg)twist-A for fixed xb- 

The twist-6 calculation indicates that the twist expan- 
sion by taking into account the multiple scattering in 
nuclear medium is actually the expansion with the pa- 
rameter asR\/Q'^ due to LPM effect. Whereas the LPM 
effect can be neglected, the expanding parameter will be 
UsRa/Q^, so that we can see the LPM effect plays an 
important role to give an additional nuclear size depen- 
dence for the parton energy loss. Furthermore, we can 
see, as is large and Ra is not very large, the param- 
eter asR\/Q^ gives a small value and it is reasonable to 
consider only the double scattering processes because all 
high twist contribution will be suppressed by the power 
oi asR\/ ■ However, for the intermediate and large 
Ra the high twist processes will give considerable contri- 
bution and the resummation of all high twist contribu- 
tion will be needed. This resummation would be similar 
to resuming all power corrections to the nuclear struc- 
ture functions in e-A DIS [23], except that in our case 
the calculations may be complicated by the additional 
induced gluon radiation. We expect the resummation of 
all high twist contributions to parton energy loss may 
adopt a similar reaction operator formalism as derived 
by GLV to compute the radiative energy loss to all or- 
ders in opacity in a hot QCD medium [3, and this work 
is in progress. 

In summary, we have calculated the semi-inclusive 
hadronic tensor and the corresponding parton energy loss 
at twist-6 in e-A DIS with helicity amplitude approxima- 
tion by utilizing the twist expansion approach. The en- 
ergy loss induced by gluon radiation at triple scattering 
(twist-6) processes has a asR\/Q^ suppression due to 
LPM interference effect as compared to the contribution 
of double scattering (twist-4) processes. The expanding 
parameter of high twist-expansion is discussed and we 
find the resummation of all high twist processes should 
be needed for intermediate and heavy nucleus with 
large nucleus radius. 
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